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ABSTRACT 



We have analyzed a broad-band optical and near-infrared spectrum of FU Tau A, a presumed young brown dwarf in the Taurus star 
forming region that has intrigued both theorists and observers by its overluminosity in the HR diagram with respect to standard pre- 
main sequence evolutionary models. FU Tau A is brighter than any other Taurus member of same or similar spectral type, and various 
phenomena (accretion, activity, binarity) have been put forth as a possible explanation. The new data, obtained with the X-Shooter 
spectrograph at the Very Large Telescope, include an unprecedented wealth of information on stellar parameters and simultaneously 
observed accretion and outflow indicators for FUTau A. We present the first measurements of gravity (log g — 3.5 ± 0.5), radial 
velocity (RV = 22.5 ± 2.9km/s), rotational velocity (vsini = 20 ± 5km/s) and lithium equivalent width (Wu = 430 ± 20 mA) 
for FUTau A. From the rotational velocity and the published period we infer a disk inclination of i >50°. The lithium content 
is much lower than theoretically expected for such a young very low mass object, adding another puzzling feature to this object's 
properties. We determine the mass accretion rate of FUTau A from comparison of the luminosities of 24 emission lines to empirical 
calibrations from the literature and find a mean of logM acc fMs/yr] = —9.9 with standard deviation a = 0.2. The accretion rate 
determined independently from modeling of the excess emission in the Balmer and Paschen continua is consistent with this value. 
The corresponding accretion luminosity is too small to make a significant contribution to the bolometric luminosity. Strong magnetic 
activity affecting the stellar parameters or binarity of FU Tau A, both combined with extreme youth, may be responsible for its position 
in the HR diagram. The existence of an outflow in FUTau A is demonstrated through the first detection of forbidden emission lines 
from which we obtain an estimate for the mass loss rate, logM out [Me/yr] < —10.4. The mass outflow and inflow rates can be 
combined to yield M out / M acc ~ 0.3, a value that is in agreement with jet launching models. 

Key words, stars: Brown Dwarfs - stars: pre-main sequence - accretion - stars: individual: FUTau A 



1. Introduction 



FU Tau A was discovered bv iLuhman e"tai1 d2009ah as the pri- 
mary component of a young brown dwarf (BD) binary in the 
Barnard215 cloud. This dark cloud is considered to be part of 
the Taurus mol ecular cloud com plex, and a distance of 140 pc is 
assigned (e.g. lTorre s et al. 2009). According to mid-infrared im- 
ages from Spitzer, only one other young star, F T Tau, is present 
within 0.5° of FUTau dLuhmanet al.ll2009bl) . Its isolated po- 
sition makes FU Tau a benchmark object for studies of brown 
dwarf formation scenarios as most mechanisms, such as e.g. 
ejection from a proto-stellar cluster, photo-evaporation, or disk 
fragmentation require th e presence of higher-mass stars, see e.g. 
IWhitworth etail (12007b. 

Luhmanetal. (2009a) obtained optical and near-IR low- 
resolution spectroscopy for the primary (FU Tau A) and an op- 
tical spectrum for the secondary. They derived spectral types 
of M7.25 and M9.25 for the two components, respectively. By 
comparison to spectral templates, they extracted differing val- 
ues for the extinction from the optical (Ay ~ 2 mag) and the 
near-IR (Ay < 1 mag) spectra of the primary, FU Tau A, and 
adopted the higher value for the calculation of its bolometric lu- 
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minosity. Masses of 0.05 M Q and 0.015 Mq were inferred com- 
paring the position of the two objects in th e HR diagram to the 
evolu tionary pre-main seq uence models of iBaraffe et al.l (1 19981) 
and IChabrier et al.l (12000b . However, as noted by ILuhman et al.l 
(2009a), in the HR dia gram FU Tau A is located well above the 
youngest (lMyr) isochrone of these models. The secondary, 
FUTauB, also seems younger than lMyr albeit by a less 
amount, contesting the general notion o f coevality for the com- 
ponent in binaries (ILuhman et alJl201QI) . 

In this article, we aim at investigating the nature of the pri- 
mary, FUTau A. Various signatures of youth are apparent in 
the available observations of this object. The presence of a cir- 
cumstellar disk was established from an analysis of its spec- 
tral energy distribution (SED), where excess emission over a 
photosphere of the same spectral type is observed in all four 
SpitzerflRAC bands, i.e. at 3.6^m, 4.5^m, 5.8/zm and 8.0^m 
(ILuhman et alJl20~09ah . Similarly, an excess in the blue part of 
the SED indicates ongoing accretion. This is bolstered by the 
high equivalent width of Ha emission measured in Luhman et 
al.'s low-resolution spectra (93 A), as well as the large width of 
the Ha pro file in a medium-res olution spectrum from Gemini 
analysed bv lStelzer etaH (12010b . The full-width at 10 % of the 
peak height was measured to be 350 km/s, distinctly higher than 
the canonical limit of 200 km/s considered to represent the bor- 
derline between magnetic activity and accretion-dominated Ha 
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emission in BDs dJavawardhana et all 20031) . Applying the cal- 
ibration provided bv iNatta et al.1 (120041) to the H a 10% width 
yielde d a mass accretion rate of 3.5 • 10 -10 M ? /vr. fStelzer et all 
(2010) caution that the Ha line is only marginally resolved in 
the Gemini spectrum. However, the mass accretion rate derived 
from the line flux of th e He 1 5876 line in the s ame spectrum, 
using the calibration of iHerczeg & Hillenbrand! (120081) is about 
a factor two higher than the value obtained from Ha, providing 
further evidence for strong accretion. 

Recently, the FUTau binary was the target of a Chandra 
X-ray observation with the aim to study magnetic activity in 
two coeval BDs of slightly different effective temperature and/or 
mass. Unexpectedly, while the secondary was not detected, the 
primary showed very strong and soft X-ray emission, reminis- 
cent of the TTauri star TWHya where the bulk of X-rays is 
produced in accreti on shocks rather than the stellar corona (e.g. 
iKastner et al.l2002h . Considering the untypically low X-ray tem- 
perature and this analogy with TW Hya, FU Tau A may be the 
first BD where X-ray emiss ion from accretion shocks has been 
detected (IStelzer et al.ll2010t) . 

The current observational picture of FU Tau presents a num- 
ber of ambiguities. First, the over-luminosity in the HR diagram, 
especially fo r the primary, allows for different i n terpre tations as 
discussed bv lStelzer etall (1201 Oh ; IScholz et al] (|2012|) . It could 
be due to extreme youth, to the primary being an unresolved bi- 
nary, or to a strong contribution from accretion to the luminosity 
of FU Tau A. Alternatively, the inhibiting influence of magnetic 
fields and/or rotation onto convection might for given luminos- 
ity lead to smaller eff ective temperature a nd, consequently, a 
wrong mass estimate ( Chab rier et al.l 12007). Indeed, the mod- 
ulations in photometric time series reveal the presence of both 
hot spots, i.e. accr etion, and cool spots, i.e. magnetic activity 
dScholz et al.ll2012T) . Secondly, the velocities derived in the Ha 
profile are smaller than the infall speeds suggested by the ob- 
served X-ray temperature, leaving a doubt on the interpretation 
of the origin of the X-ray emission in accretion shocks. 

Aiming at a better understanding of this puzzling BD, we 
have obtained broad-band spectroscopy from the UV to the near- 
IR and multi-color time-series photometry for FU Tau A. The 
motivation for collecting these data was to seek for a better un- 
derstanding of its accretion and activity characteristics as possi- 
ble cause for the overluminosity of FU Tau A in the HR diagram 
through the study of spectral signatures. The new photometry is 
part of our study of FU Tau A's long-term variability. The obser- 
vations and data analysis are described in Sect. [2] Stellar prop- 
erties are derived in Sect. [3] In Sects. |4]and|5]we examine the 
outflow and accretion characteristics of FU Tau A, and in Sect. [6] 
we discuss our results. 

2. Observations and data reduction 

FUTau A was observed on Jan 11, 2011 with the X-Shooter 
spectrograph at t he VLT. The data were acquired within the 
INAF/GTO time dAlcala et alj201 ll) . With its three spectrograph 
arms, X-Shooter provides simultaneous wavelength coverage 
from 300 - 2480 nm. Slit widths of 1.0"/0.9"/0.9" were used 
in the UVB/VIS/NIR arms, respectively, yielding spectral reso- 
lutions of 5100/8800/5600. The total exposure time in each of 
the three spectrograph arms was 1800 sec. We obtained a signal- 
to-noise ratio of 1 - 8 in the UVB, 10 - 20 in the VIS, and 20 - 30 
in the NIR arm. The data were obtained in nod mode and were 
reduce d independently for eac h arm with the X-Shooter pipeline, 
vl.3.7 (Modig hani et all2010l) . Following the standard steps in- 
cluding bias or dark subtraction, flat fielding, optimal extrac- 



tion, wavelength calibration, sky subtraction, correction for at- 
mospheric extinction and flux calibration. However, the pipeline 
flux-calibrated spectra are not corrected for telluric absorption 
bands. The telluric correction was performed independently in 
the VIS and NIR spectra in the following way. For the VIS arm 
the spectrum of the telluric standard (HIP 20789, spectral type 
B7 VY, normalized to the continuum, was used as input in the 
IRAFQ task "telluric". For the NIR spectrum, a response func- 
tion was first derived by dividing the non flux-calibrated spec- 
trum of the telluric by a black-body of the same effective tem- 
perature (T o ff=13,000K) as the telluric standard. This response 
function, containing the telluric lines, was then used as input in 
the IRAF task "telluric". In this way, the telluric correction and 
the correction for the response function were done simultane- 
ously. Although the shape of the resulting spectrum after this 
procedure is correct, the flux calibration is only relative to the 
response function. To bring the NIR spectrum into an absolute 
flux scale we multiplied with a factor that was estimated using 
the pipeline flux-calibrated spectrum. Finally, wavelength shifts 
due to instrumental flexures were corrected using the flexcomp 
package within the pipeline. The precision of the wavelength 
calibration is better than 0.01 pix corresponding to 0.002 nm in 
the UVB and VIS arms and to 0.006 nm in the NIR ami. 

The flux-calibrated spectrum of FUTau A for the full X- 
Shooter wavelength range is shown in Fig. Q] Photometry from 
the literature is overlaid, demonstrating the high quality of the 
flux calibration. In particular, the match of the different instru- 
mental arms is very good. In fact, during the observation the 
seeing was excellent (« 0.85"). 

3. Stellar properties 

3. 1. Spectral type and extinction 

We estimated the spectral type and extinction of FU Tau A by 
comparing it to the spectral templates defined by Manara et al. 
(2013, A&A subm.) on the basis of X-Shooter spectra for 24 
non-accreting and unabsorbed YSOs in various star forming re- 
gions. These spectra define a continuous spectral sequence from 
M0 to M6.5 with steps of 0.5 in spectral subclass. In addition 
they comprise two objects at the end of the M subclass. At their 
young ages (1 — lOMyr) these objects are expected to have simi- 
lar gravity to FU Tau A. We determined simultaneously the spec- 
tral type and the extinction of FU Tau A by artificially redden- 
ing the templates between Ay = 0...2 mag until the best match 
to FU Tau A was found. For the reddenin g we used the extinc- 
tion law of Weing artner & Drainel (12001 1) . In Fig. [2] the spec- 
trum of FU Tau A is compared to that of Par-Lup3-1 (M6.5) and 
DENIS-PJ124514.1-442907 (M9). The spectrum of Par-Lup3- 
1 appears very similar to FUTau A and we conclude that the 
spectral type of FUTau A is M6.5 or slightly later. The gap in 
the spectral type sequence of Manara et al. (2013) does not al- 
low us to put stronger constraints. We estimated an extinction of 
Ay ~ 0.5 ± 0.5 mag. The templates in Fig.|2]are both reddened 
by this amount. 

In a second approach to determine the spectral type of 
FU Tau A we have c alculated various spectral indices following 
Ridd ick et alj d2007l) for optical w avelengths, and the H 2 in- 
dex defined bvlAllers et alj y007) and the H2O-K2 index from 
Roias -Avala et al.l (120121) for the near-IR spectrum. The indices 
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Fig. 1. Broad-band flux-calibrated X-Shooter spectrum of FUTau A. Published photometry from iLuhman et aD (l2009al) . from 
IScholz et al.l d2012l) and unpublished RIJK photometry obtained in September 2012 with Andicam at the 1 .3m telescope on Cerro 
Tololo as part of the SMARTS program DUBLIN- 1 2 A-00 1 is overplotted. 
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Fig. 2. X-Shooter VIS spectrum of FUTau A compared to two 
spectral templates from Manara et al. (2013). A reddening of 
Ay = 0.5 mag has been applied to the templates. The spectra 
are normalized at 750 nm and the templates are vertically offset. 



The spectral type we obtained from the X-Shooter spectrum 
with the two methods explained above is also consistent with 
our earlier results from low-resolution spectroscopy where we 
found values betwe en M6.6 and M6.8 in a series of five spectra 
(Scholz et al. 2012). Using a different set of spectral templates 
ILuhman et alJ (l2009ah found a similar, slightly cooler spectral 
type ofM7.25± 0.25. 

Luhman et al.'s value for the extinction derived from their 
optical spectrum (Ay < 1 mag) is also similar to ours but they 
found a higher value from the near-IR spectrum (Ay = 2 mag). 
The latter value is clearly incompatible with the X-Shooter NIR 
spectrum for which we derive Ay ~ 0.75 ± 0.5 mag in an anal- 
ogous way as the one described for the VIS. Throughout this 
paper we adopt the spectral type and extinction derived from 
the VIS spectrum. We expect the determination of extinction us- 
ing zero-extinction spectral templates observed with the same 
instrument to be preciser in the VIS than in the NIR because, 
generally, the effects of extinction are larger in the VIS than in 
the NIR. 



that we have used are consistent with the spectral sequence of 
Manara et al. (2013). The flux ratios for the spectral indices 
were calculated from the spectrum of FU Tau A after modifying 
it by various amounts of extinction between Ay = and 2 mag. 
We co nfirm that the optical spectral indices from iRiddick et alJ 
(2007) are independent on extinction for the range of extinc- 
tions tested here and derive a spectral type of M6.6 ± 0.4 for 
FU Tau building the average and standard deviation of the V02, 
Rl, R2, R3, Ti08465, and V07445 indices. The H 2 index 
yields M6.7 ± 0. 2 and the H 2 Q-K2 in dex M7.1 for a range of 
Ay = 0...1 mag. IRiddick et all d2007h and lAllers et al.l d2007l) 
provide calibrations specifically for young M dwarfs, i.e. their 
spectral indices are independent on gravity, and we find excel- 
lent correspondence in the derived spectral types. The slightly 
later spectral type o btained from the H9O-K2 index may be re- 
lated to the fact that lRoias-Avala et alJ d2012l) have calibrated it 
for nearby (evolved) M dwarfs. 



3.2. Stellar parameters 

Our new measures for the spectral type and the extinction can be 
used to revise the stellar parameters of FUTau A. The radius of 
FUTau A can be obtained from Stefan-Boltzmann's law. Rather 
than using the literature value for the bolometric luminosity we 
recompute it taking account of the difference i n our a dopted 



Ay and the value assumed by Luhman etal. ( 2009al) . From 
the J magnitude given by Luh man etal. ( 20*09al) and a bolo- 
metric correction of B.C.i = 2.05 dDahn et alJl2002l) we find 
^boi = 0-13 Lp). We have m easured a slightly earlier spectral 
type than ILuhman et al] d2009ah . yielding an effective tempera- 
ture of T c ff = 2940 K fo r a temperature scale intermediate be- 
tween dwarfs and giants dLuhmanlll999t) . For these numbers we 
derive i?* = 1.4 i? Q , substantially lower than th e value result- 
ing fr om the luminosity and temperature given by Luhma net al.l 
d2009al) . i?* = 1.8 R G . Extrapolating the position of FUTau A 
in the HR diagram down to the youngest isochrone of the evolu- 
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Fig. 3. Detail of the spectrum of FU Tau (black solid lines) in 
the wavelength range around the Kl (upper panel) and Nal 
(lower panel) absorption doublets. The spectrum of FUTau has 
been corrected for telluric absorption lines. The red dotted line 
and the green dashed line represent synthetic spectra with log g 
as labelled for T e g = 2900 K and rotationally broadened to 
usini = 20km/s. 



tionary models bv IChabrier et al.l (120001) . the mass is 0.08 M , 
placing FU Tau A right at the hydrogen burning mass limit. 

We used the gravity- and temperature-sensitive absorption 
doublets of Nal at AA 818.33, 819.48nm and of Kl at AA 
766.48, 769.89 nm to determine logg by comparison of the 
X-Shooter data to synthetic spectra. We retrieved BT-DUSTY 
model spectra for a range of T e g around the expected value of 
2940 K and a range of log g values from the star, brow n dwarf & 
planet atmosphere web simulator dAllard et al.ll20 10). electroni- 
cally availably In Fig.[3]the synthetic spectra for T eS = 2900 K 
with three different values for log g are overlaid on the spectrum 
of FU Tau A in the region of the Na and K doublets. The model 
spectra have been rotationally broadened to vsini = 20km/s, 
the rotation rate determined for FU Tau A in Sect. 13.51 The width 
of the observed lines in Fig.[3]is in good agreement with a gravity 
of log g = 3.5 ± 0.5. The uncertainty is due to the 0.5 steps for 
log g in the grid of synthetic spectra. The gravity derived from 
the spectrum a grees with the value expected from the evolution- 
ary models of Bara ffe et al.l ([1998) adopting the stellar parame- 
ters derived above and an age of 1 Myr. 

3.3. Lithium absorption 

Theoretical models predict that the lithium content of low-mass 
stars gets rapidly depleted throug hout the first ~ 50 Myrs of the 
pre-main sequence evolution (e.g. lD'Antona & M azzitelli 1997; 
Baraffe et al. 1998). In particular, the LiIA670.8nm absorption 
line is a well-known age indicator in young low-mass objects. 
This was shown, e.g., in studies of the line equivalent width 
for c lusters and associa tions with different HR diagram ages 
(e.g. iMentuchet all 120081) . The age at which lithium depletes 
increases with decreasing mass. For the fully convective very 
low-mass (< 0.35 Mq) objects in star forming regions standard 
evolutionary models predict that the original abundances a re still 
retained until an age of at least 10 Myrs (e.g. Jeffries 2006) and a 
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Fig. 4. Li lA670.8nm absorption in FUTau A (top) and TiO ab- 
sorption band at 615 nm, a diagnostic of veiling, compared to 
the same spectral region in the unveiled Class III template Par- 
Lup3-1 (bottom). 



cosmic abundance of A(Li) « 3.1 (Zap atero Osorio et al . 2002) 
is expected. 

The Li I A670.8nm absorption of FUTau A is shown in 
the top panel of Fig. [4] We determined an equivalent width of 
0.43 ± 0.02 A. In practice, for accreting objects photospheric 
absorption lines may be filled in by excess continuum emis- 
sion leading to shallower lines and an underestimate of the line 
equivalent widths. The amount of this 'veiling' can be estimated 
by comparison with the spectrum of an unveiled star of the 
same spectral type. The only available template for FUTau A 
observed with the same instrument and the same spectral reso- 
lution is Par-Lup3-1 (M6.5). We have thus measured the EWs 
of the T iO absorption band s in the region between 660 nm and 
680 nm dValenti etalj 19981) for both the template and FU Tau A. 
We calculated the ratio -EWtempiate/ EWfut^uA for eight dif- 
ferent bands and found that the ratio is at most 1.2. A more 
detailed evaluation is prohibited by the noise level of the tem- 
plate spectrum. If these differences in EWs are due to veiling 
the corrected EW(Li) value for FU Tau A would be < 516 mA. 
We have repeated the same exercise using the spectral features 
shown in the bottom panel of Fig. [4] The result is a correction 
factor of < 1.3 which is consistent with weak veiling. 

In Fig. [5] we compare our measurement to other studies of 
Uthium in YSOs of late-M sp ectral t ype from the litera ture: <j Ori 
( Zap atero Osorio et al.l2002l) . A Ori dBavo et al.l201 ll). and three 
very low-mass (VLM) objects in Taurus (Barrado v Navascu6sl 
2004). A compilation of lithium equ ivalen t width measurements 
in Taurus is given bv ISestito et al.l d2008l) for a larger range of 
spectral types. We do not consider this sample here because it in- 
cludes a number of stars with dubious membership. For the case 
of A Ori we consider only the subsample with confirmed clus- 
ter membership and with spectroscopically determined spectral 
types. 

Accreting objects are highlighted with filled plotting sym- 
bols in Fig. [5] and have been selected as follows: In cOri we 
consider accreting all objects with Wh„ > 20 A, in A Ori we 
rely on the 'accretor flag' given by Bav o et al.l (1201 lb which 
is based on Ha emission and in Taurus we consider accret- 
ing objects with W^ a w% > 300km/s. The choice of the 10% 
width rather than the equivalent width as accretion diagnostic for 
Taurus is motivated by the fact that one object, KPNO-Tau 5, has 
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Fig. 5. Lithium equivalent width versus spectral type for 
FUTau A (filled star symbol) and other measurements of M- 
type YSOs collected from the literature (see text in Sect J3.3l >. 
Accretors are characterized by filled symbols. Also plotted are 
the curves of growth for log g = 4.0 and three different abun - 
dances with labels, calculated by Zapa tero Osorio et al . (2002). 



Who = 21.1 A, at the borderline to the accretion regime but is 
classified as non-accretor due to the narrow Ha profile and weak 
He I A 667.8 emission. 

The above-mentioned star forming regions all have an age of 
1 — 5 Myr and the lithium content of M stars is expected to re- 
flect the initial cosmic abundance, A(Li) = 3.1. However, a large 
spread of equivalent widths is observed for a given spectral type. 
Observations obtained with different spectral resolution may re- 
sult in systematic differences of the line measurements but these 
errors are smaller than the observed spread of the equivalent 
width values. This spread also corresponds to an unexplained 
spread in abundances as evident by the over-plotted curves of 
growth for log g = 4.0 calculated by IZapatero Osorio et al.l 
(2002). FUTau A presents the lowest Li measurement observed 
so far for its spectral class. Taking into account veiling, FU Tau A 
moves up in the diagram but is still well below the theoretically 
expected abundance. 

3.4. Radial velocity and UVW velocity 

We have used the strongest absorption lines, the KI doublet 
at AA 766.48, 769.89nm and the Nal doublet at AA 818.33, 
819.48 nm, to estimate the radial velocity of FU Tau A by com- 
parison of the observed line centers to the rest wavelengths 
extracted from the NIST Atomic Spectra Database^ After ap- 
plication of the barycentric correction (— 20.6km/s) the mean 
Doppler shift of these lines yields RV = 22.5 ± 2.9km/s 
where we have added to the standard deviation of the four 
measured lines the uncertainty of the wavelength calibration 
(see Sect. [2). The local-standard-of-rest velocity of FU Tau A is 

12. 5±2.9km/s. 

iBertout & Genov a (2006) have compiled a list of 127 stars in 
the Taurus star forming complex with RV measurements. These 
authors note that no systematic high-precision RV survey has 
been performed in Taurus. Their RV distribution is strongly 
peaked at 15 — 16km/ s, in agreement w ith histor ic smaller 
samples presented by Hartmann et al. (1986|) and IWalter et al.l 

3 http://www.nist.gov/pml/data/asd.cfm/ 



Table 1. Proper motion of FUTau A from the literature and 
galactic velocity calculated using our RV measurement and 
rf=140pc. 



Reference 


Ma 


,fJ.s 


{U, V, W) 






[mas/yr] 


[km/s] 




Ducourant et al. (20051 


14.7,- 


-26.0 


(-22.09,-9.08,- 


-10.36) 


Luhman et al. (2009a) 


7.2,- 


-17.5 


(-23.54,-16.55,- 


-10.41) 
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Fig.6. The Nal absorption doublet AA 818.33, 819.48 nm of 
FU Tau A is shown as continuous black line. The synthetic spec- 
trum, broadened at projected rotational velocities of 15, 20 and 
30 km/s, is represented with the green-dotted, red-dashed and 
blue-dotted lines, respectively. 



(1988). About ~ 5 7r of this list have R\ ' > 22 km/s, i.e. the R\ ' 
of FU Tau A is only marginally compatible with that of Taurus. 

We have combined our RV measurement with the proper 
motion of FUTau A given in the literature to obtain the 
UVW velocities. We use both th e proper motion valu es of 
Ducourant et al. 1 (120051) and those o f iLuhman et all (l2009al) . The 
resulting galactic velocity components for a distance of 140 pc 
are listed in Table Q] The average space motion for Taurus 
for the same distance is (U,V,W) = (-16.5 ± 4.6, -13.2 ± 
2.5,-11.0 ± 4.0) km/s ( IBertout & Genovall2006l) . FUTau A is 
in rough agreement with these values. 

3.5. Rotational velocity 

The projected rotational velocity, vs'mi, was estimated by 
comparing the profile of the Nal absorption doublet at 
AA818.33,819.48nm with that of a synthetic BT-DUSTY spec- 
trum of the same effective temperature and gravity as FU Tau A. 
The synthetic spectrum was gathered f rom the star, brown dwarf 
& planet atmosphere web simulator (All ard et al.ll2010l) using 
a v sini — Okm/s. In order to reproduce the profile of the Na I 
absorption doublet in FU Tau A, the synthetic spectrum was con- 
volved with rotational profiles (see Gray 1992) of several values 
of v sin i. The best match is for usini = 20 km/s (see Fig. [6j>. 
We estimate an error on the order 5 km/s. 

We can combine our measurements of the rotational veloc- 
ity ( v sin i) with the phot ometrically determined rotation period 
from lScholzetalJ d2012l) (P = 4.0 ± 0.2 d) and find Rsini = 
1 .6 ± 0.5 R©. The value obtained for R sin i is compatible with 
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our new value for the stellar radius = 1.4 Kq) within the 
errors and implies a disk inclination angle i > 53°. 

4. Outflows 

We have searched for forbidden emission lines (FELs) as sig- 
natures of shocks forming in outflows. A list of the equivalent 
widths (EWs) and line fluxes for the FELs detected in FU Tau A 
is given in Cols. 3 and 4 of Table [2] The EWs and fluxes are 
obtained from an by-eye estimate of the local continuum. The 
uncertainties of the EWs and the fluxes represent the mean and 
standard deviation from three measurements carried out on each 
line. In our approach of estimating these uncertainties we take 
account of the fact that the uncertainties in the flux are domi- 
nated by the uncertainty of the extinction while the major source 
of the uncertainties in the EWs are the statistical fluctuations. 
Therefore, for the errors of the EWs we have measured the line 
in the spectrum dereddened by 0.5 mag at another two positions 
corresponding to our estimate of the upper and lower value of 
the adjacent continuum emission and computed the standard de- 
viation of the three measurements. For the uncertainties of the 
line fluxes we have measured the line in the observed spectrum 
(corresponding to an assumption of Ay = Omag) and in the 
spectrum dereddened for Ay = 1 mag, and combine those two 
values with that from the spectrum dereddened by 0.5 mag to get 
the standard deviation. 

The non-detection of th e rSIIlA671.6nm line does n ot allow 
us to apply the technique of lBacciotti & Eisi offel ( 1999) for esti- 
mating plasma parameters from the line ratios of FELs. For a de- 
tailed analysis of luminosities, kinematics and mass outflow rate 
(M out ) of FU Tau A we consider only the two most prominent 
lines, [OI] A630.0 nm a nd fSIIl A673. 1 nm, an d we follow the pre- 
scription described by Harti gaiTet al.l (1 1995b who presented re- 
lations of the type 

logM out = ^■(l+-)( T F7T^)(^^)- 1 ( : r^)(l) 
n e 150 kms 1 2-10 i!5 cm L sun 

The numerical constant X is derived in Hartigan et al. (1995) 
for the [SII]A673.1nm and the [OI]A630.0nm lines. Eq.HJin- 
volves the electron density (n c ), the tangential velocity of the 
outflow (Vt) and the projected size of the aperture in the plane 
of the sky (It). For a 1" slit we obtain for the distance of Taurus 
It = 2 ■ 10 15 cm. For the reasonable assumption that the outflow 
is perpendicular to the disk we use the minimum possible disk 
inclination angle (idisk ~ 50°; see Sect. 13.5b and the measured 
RV of the FELs with respect to the stellar RV to estimate the 
tangential velocities. The derived quantities are given in Table[2] 
The uncertainty in the velocities are given in col. 7. They include 
the standard deviation of three wavelength measurements and 
the uncertainty of the stellar RV. 

We have no measurement of the electron density for the out- 
flow of FU Tau A. A lower limit to n c is given by the fact that 
we detect [SII] A673.1 nm but not [SII]A617.7nm emission. The 
flux ratio between these t wo lines is, theref ore, larger than unity 
and n c > 10 3 cm~ 3 (e.g. lOsterbro ck 1989). This gives an upper 
limit to the mass loss rate determined from [SII] A673. 1 nm of 
log M Q \xt [Mo /yr] < — 10.4. Assuming the same lower limit for 
the electron density, the [OI]A630.0nm line gives a much less 
stringent constraint on Mout- We note that the mass loss rates 
computed from the two lines agree for n c ~ 10 5 cm -3 . For that 
case, we find logM out [M /yr] < —11.5. Recall also that the 
value of the inclination angle assumed in this calculation is the 



minimum possible value on the basis of the rotational proper- 
ties of FU Tau A. If the disk inclination is higher than that value, 
the upper limit on the mass loss rate is higher and less meaning- 
ful, e.g. for i — 75° the mass loss rate obtained from [SII] is 

< 10" 10 M o /yr for n c > 10 3 cm- 3 and< 10" 110 M Q /yr 
for n e ~ 10 5 cm~ 3 . 



5. Accretion 

Different methods have been described in the literature for es- 
timating the mass accretion rate (M acc ) from the properties 
of emission lines or continuum emission. The various accre- 
tion diagnostics probe several physical regions, e.g. the ex- 
cess continuum that veils the absorption fe atures of classical T 
Tauri stars is ascribed to accretion shocks (ICalvet & Gullbrin^ 
1998) and the high-excitati on Hel emission to post-shock re- 
gions dBeristain et al J 120011) . Hydrogen lines from the Paschen 
and Balmer series are produced in accretion flows, with Ha, hav- 
ing the largest optical depth, originating in th eir outermost parts 
(iHartmann et al.lll994tlMuzerolle et alfl998l) . 

Here, we use two approaches for measuring the mass accre- 
tion rate of FU Tau A: We make use of empirical relations with 
emission line fluxes and luminosities (Sect. lBTXH and we model 
the continuum excess (Sect. l5T2~2l ). 

5. 1. Emission line analysis 

Using IRAF we have measured equivalent widths and line fluxes 
of the dominant emission lines in the spectrum in the same man- 
ner as described in Sect.|4]fbr the case of the FELs. 

The by far strongest line in the X-Shooter spectrum of 
FU Tau A is Ha. We measure an equivalent width of 93 ± 5 A. 
Previous Ha measurements of FU Tau A reported in the litera- 
ture, all b ased on lower-resolut i on spectra, range between 93 and 
155A (iLuhman et alJ l2009bt Isteizer et al.l l2010t IScholz et all 
12011 . Considering that low-resolution observations tend to un- 
derestimate line equivalent widths, this points at only modest 
variability. The Balmer series can be reliably identified up to 
n = 16. No lines of the Paschen or Brackett series are detected. 

Some of the Balmer line profiles are displayed in Fig. [7] All 
Balmer lines including Ha are single-peaked. The vertical lines 
in Fig.|7]represent the expected line center corrected for barycen- 
ter motion and the systemic RV of FU Tau A derived in Sect. 13.41 
While the Ha profile is nearly symmetric the other low-n Balmer 
lines show small asymmetries with a flux deficiency on the red 
side of the profile. This mi ght be interpreted as r ed absorption, a 
typical signature of infall (Hartman n et al.|[l994l) . 

5.2. Mass accretion rates 
5.2.1. Emission line diagnostics 

The published calibrations between emission line strength and 
mass accretion rate come in two flavors. The first consists in re- 
lations between line luminosity (Xii„ c ) and accretion luminosity 
(L acc ) of the type 

log (L acc /L Q ) = a ■ log (Liino/L©) + b (2) 
This can be converted to accretion rate according to 

M acc = 1.25-^-^ (3) 
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Table 2. Forbidden emission lines: nominal wavelength (col. 2), line equivalent width and flux (cols. 3 and 4), radial velocity with 
respect to the barycenter-corrected stellar velocity (col. 5), tangential velocity (col. 6) and uncertainty of the velocities. 



Element \ EW / line • 10 iB KV Vt Al7 





[nm] 


[nm] 


[erg/cm 2 /s] 


[km/s] 


[km/s] 


[km/s] 


SII 


406.860 


-0.305 ±0.102 


2.7 ± 1.5 


-26.6 


99.2 


3.2 


SI1 


407.635 


-0.081 ± 0.083 


0.4 ±0.2 








Ol 


630.030 


-0.323 ±0.053 


8.8 ±3.4 


-12.6 


47.2 


3.2 


OI 


636.378 


-0.067 ±0.011 


2.7 ± 1.0 








Nil 


658.260 


-0.009 ± 0.007 


0.7 ±0.4 








SII 


673.082 


-0.030 ±0.010 


1.2 ±0.4 


-11.0 


41.1 


3.2 






Fig. 7. Line profiles for the lower Balmer series of FU Tau A. The dashed vertical line represents the expected line center correcting 
for the barycenter motion and the stellar RV computed from absorption lines. 



where the numerical factor deri ves from the ass umption for the 
inner disk radius, i? in = 5i?* dGullbring et al.| [T998). The sec- 
ond approach relates the line surface flux (inline) directly to the 
accretion rate in the form of 

log M acc = C ■ log inline + d (4) 

The coefficients a, b, c and d have been derived for an increasing 
sample of YSOs, however, so far without consensus between the 
various studies. 

To estimate M acc for FU Tau A we have used t h e pub - 
lished relations pre sented by iHerczeg & Hillenbrand! (|2008), 
iMohantv etal] d2005l) . Rigliaco et al. (2012, A&A in press; here- 
after RNT12) as well as new calibrations by Alcala et al.,(in 
prep). The latter ones are based on an X-Shooter sample of 
36 YSOs from the Lupus star forming region and represent the 
largest homogeneous analysis of this kind performed so far. The 
distance used to convert line fluxes to luminosities is that of 
the Taurus star forming complex, 140 pc. The mass and radius 
of FUTau A required for evaluating Eqs. [3] and |4] are given in 



Sect. 13.21 The derived mass accretion rates for all detected emis- 
sion lines in common with the relations presented in the litera- 
ture are shown in Fig. [8] Each literature source is marked with 
a different plotting symbol. Some values that deviate strongly 
from the bulk of the measurements are not shown for clarity but 
discussed below. The error bars represent again the uncertainty 
in the reddening with upper and lower bounds corresponding to 
Ay = and Ay = 1 mag, respectively. 

We have calculated the mass accretion rate for all emission 
lines detected in the X-Shooter spectrum of FU Tau A that have 
been calibrated in the literature. However, we have computed 
the average of M acc only from those lines that we consider the 
most reliable accretion diagnostics. In particular, we excluded 
the Ca II IRT lines and Na I D, for which strongly discrepant ac- 
cretion rates are found with different calibrations. Some of the 
values for M a cc that we have derived for these lines are outside 
the range of values shown in Fig. [8] For a detailed investiga- 
tion of the reliability of individual emission lines as accretion 
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diagnostic see e.g. Alcala et al. (in prep.). We also did not take 
into account [Ol]A630nm and Ha that are often affected by 
winds in YSOs (see Sect. 0J. CallK and He are also excluded 
because these lines are partially blended in our spectrum. The 
average obtained from the remaining lines shown in Fig. [8] is 
(log M^) [M /yr] = -9.9 ± 0.2 for 24 lines from Alcala et 
al., (logM acc )[M /yr] = -9.8 ± 0.2 for 6 lines from RNT12, 
and <logM acc )[M /yr] = -10.0 ±0.3 for 11 lines from HH08. 
Here the uncertainties represent the standard deviation of all 
measurements. The mean and standard deviation of the calibra- 
tions by Alcala et al. are overplotted in Fig. [8] The influence of 
the uncertainties of mass and distance on this result is examined 
in Sect. [6] We anticipate here that they make log M acc change 
only marginally. Similarly, a smaller value for the inner disk ra- 
dius, e.g. i?i n = 2 Rq, would yield a mass accretion rate that is 
marginally compatible with the value given above. 

We have also derived the mass accretion rate from the Ha 
10 % width. The measured full width half maximum at 10 % 
of the peak heig ht is 228km/s. For the calibration provided by 
iNatta et al] d2004l) this corresponds to a mass accretion rate of 
logM acc [M /yr] = —10.7 ± 0.5. This value, also shown in 
Fig. [8] is almost one order of magnitude lower than the mean 
value derived from the relations presented by Alcala et al. 

5.2.2. Continuum excess 

The accretion luminosity can be estimated directly by measuring 
the continuum emission in excess of the photospheric one. This 
is done by comparison to a non-accreting template star of the 
same spectral type. Unfortunately, the Class III Par-Lup3-1 that 
we have used as spectral template in Sect. 13. H is too noisy to be 
used for this purpose. Fig.|9]shows the blue part of the spectrum 
of FU Tau A compared to that of another young Class III star of 
similar (albeit slightly earlier) spectral type, SO 925 (M5.5; see 
RNT12 and Manara et al. 2013) normalized in the region 700 — 
720 nm. FUTau A shows a clear excess emission in the Balmer 
continuum and a relatively small Balmer jump, defined as the 
ratio of the fluxes observed in the continuum at 360 nm and at 
420 nm, of ~ 1.2. The Balmer limit is at 346.6 nm but due to the 
blending of the higher Balmer lines it appears shifted to 370 nm. 

We model the excess emission as due to a slab of hydro- 
gen of given electr on density, temperature and length, (see e.g., 
IValenti etal]|1993l) : HH08; RNT12. and obtain an estimate of 
the excess continuum luminosity of L acc ~ 1.3 x 1O~ 4 L . 
Fig- shows the best-fitting model in blue. It reproduces quite 
well the observed Balmer continumm, the Balmer jump and the 
profile of the Ca II line at 420 nm. Veiling of photospheric lines 
at longer wavelengths is very weak as shown in Sect. 13.31 The 
uncertainties of L acc are quite large (50 % at least), due to the 
noise of the spectra of both FU Tau A and the Class III template 
and due to the mismatch between the spectral type of FU Tau A 
and that of the template. For M = 0.08 M , = 1.38 R 
as determined in Sect. I3.2l the corresponding mass accretion rate 
from Eq. [3] is logM acc [M /yr] ~ —10.1, in reasonably good 
agreement with the estimates from the emission lines presented 
in Sect. EXT] 

6. Discussion 

We have analysed a broad-band (300 — 2480 nm) medium- 
resolution (R ~ 5000 — 9000) spectrum of FU Tau A obtained 
with X-Shooter. On the basis of these data we have presented the 
first measurements of gravity, radial velocity, rotational veloc- 




flim II i I' III i ' I L I I I il I I I I I I I 

340 360 380 400 420 

Wavelength (nm) 



Fig. 9. Balmer jump region in the spectrum of FUTau A, heav- 
ily smoothed to reduce the errors on the continuum (red), the 
Class III template SO 925 (black), the excess emission from the 
accretion slab model (dashed green) and the resulting modelled 
emission, i.e. template plus excess emission, in blue. All spec- 
tra are normalized in the region 700 — 720 nm. Note that the 
observed emission in the range 346.6 nm to ~ 370 nm is domi- 
nated by unresolved Balmer lines, which are not included in the 
slab model. 



Table 3. Equivalent widths and fluxes for emission lines. 



Line 


Ao 


EW 


log /li nc 




[nm] 


[nm] 


[erg/cm 2 /s/nm] 


Hel 1083 


1083.038 


-0.150 ±0.026 


,0 7 q + U.(J« 

lO. IO_q .07 


Call 866 


866.381 


-0.017 ±0.003 


-15.11Io.l? 


Call 854 


854.347 


-0.076 ±0.015 


-14.611°;!° 


Call 849 


849.936 


-0.042 ±0.010 


-14.861™ 


Hel 667 


667.928 


-0.089 ±0.014 


-15.401™ 


Ha 


656.374 


-9.766 ± 0.709 


-13-181™ 


[OI]630 


630.093 


-0.314 ±0.039 


-15.09j|g 


NalDl 


589.697 


-0.221 ±0.044 


— 15.51_ 21 


Nal D2 


589.096 


-0.372 ± 0.067 


-15.271™ 


Hel 587 


587.670 


-0.581 ± 0.065 


-15.011™ 


Hel 501 


501.625 


-0.071 ±0.016 


-16.071™ 


H/3 


486.184 


-6.360 ± 1.765 


-14.141^ 


Hel 471 


471.397 


-0.031 ±0.009 


-16.50±g;SS 


Hell 468 


468.633 


-0.048 ±0.016 


-16.301™ 


Hel 447 


447.216 


-0.470 ±0.164 


-15.471™ 


H 7 


434.096 


-6.334 ± 2.929 


-14.441™ 


H6 


410.221 


-4.151 ± 2.327 


-14.641 ,; 2 . 2 


Hel 402 


402.687 


-0.238 ± 0.120 


-15.871™ 


He 


397.068 


-1.767 ± 1.078 


-14.891™ 


CallH 


396.898 


-5.297 ± 3.764 


-14.63l°'j« 


CallK 


393.419 


-5.981 ±2.717 


-1452lg;» 


H8 


388.944 


-4.875 ± 3.823 


-14.781^ 


H9 


383.588 


-2.392 ± 1.852 


-14.98123 


H10 


379.845 


-0.975 ± 0.447 


-15 1^+ - 24 


Hll 


377.113 


-0.855 ±0.512 


15.25 o 27 


H12 


375.071 


-1.230 ±0.684 


-15.281™ 


H13 


373.484 


-0.658 ± 0.422 


-15.44lg; a i 


H14 


372.242 


-0.384 ± 0.250 


-15.66+gi? 


H15 


371.232 


-0.377 ± 0.276 


-15-641™ 


H16 


370.432 


-0.105 ±0.057 


-16.04!°,;*? 
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Fig. 8. Mass accretion rate for FUTau A derived from various emission line fluxes and luminosities and from the 10% width of 
Ha using different calibrations provided in the literature as labeled in the figure. Dashed and dotted lines represent the mean and 
standard deviation for the relations from A13 using 24 emission lines. 



ity and lithium content as well as the first detection of outflow 
activity in FUTau A. Moreover, this spectrum provides a large 
number of accretion diagnostics which we have employed for a 
detailed evaluation of the mass accretion rate. 

The motivation for our spectroscopic study was the search 
for an explanation of the (apparent) overluminosity of FU Tau A 
with respect to the predictions of evolutionary models for its 
age and mass. This requires precise measurements of the fun- 
damental stellar parameters (T e ff, -^boi, R*, loggO as well as a 
study of other phenomena such as accretion and magnetic activ- 
ity that may lead to a wrong estimate of those parameters. We 
have discussed several possible factors that may influence the 
position of FU Tau A in the HR diagra m in our previous papers 
dStelzer et al.l l2010; Sch olz et al.ll2012l) and we resume this dis- 
cussion here, adding the wealth of information provided by the 
X-Shooter spectrum. 

6. 1. Accretion and outflow in FU Tau A 

In our detailed investigation of mass accretion in FU Tau A we 
made use of an unprecedented large number of empirical rela- 
tions between the luminosity of individual emission lines and the 
accretion luminosity derived from a comprehensive X-Shooter 
study of accretors in the Lupus clouds (Alcala et al., in prep.). 
These relations comprise emission features collected with all 
three arms of the X-Shooter spectrograph. As compared to sim- 
ilar relations presented in the literature, we consider them the 
most reliable calibrations for our observation of FU Tau A given 
that they were derived with the same instrument. A mean value 
of (logM acc )[M /yr] = -9.9 ± 0.2 is found from a total of 
24 emission lines. The accretion rate obtained from the HalO % 
width (log M acc lo o/ o [M /yr] = — 10. 7±0. 5) is lower by almost 

one order of magnitude. The discrepancy of log Af acc 10 % with 
respect to other accretion diagnostics supports previous notions 
of the poor reliability of this tracer (e.g., Costigan et al. 2012, 
MNRAS in press). On the other hand, the accretion rate from 
the Hal0% width as measured in the X-Shooter spectrum is 
also much lower than the value measured with the same diagnos- 
tics fr om a previous Gem ini spectrum (logM acc l0 % [Mq/vi'] = 
— 9.5: IStelzer et al.l2010l) . The low velocities in the new Ha pro- 
file are not sufficient to explain the temperature of the X-ray 



emitting plasma observed by Stel zer et all (1201 Ol) as originating 
in an accretion shock. 

We have presented the first evidence for outflow activity in 
FUTau A in the form of several forbidden emission lines. The 
outflow of FU Tau A has recently also been detected at millime- 
ter wavelengths (Monin et al., A&A subm.). In this paper, we 
have used the [OI] A630.0nm and [SII] A673.1nm lines de- 
tected in the X-Shooter spectrum to derive the mass loss rate 
of the outflow using the approach presented by Hartiga n et al.l 
(1995). We trust the result from [SII] A673.1nm more than 
that for [OI] A630.0nm because the former has a weaker de- 
pendence on the unknown electron density. An additional un- 
certainty is introduced by the fact that we have only a lower 
limit on the disk inclination angle that determines the transverse 
outflow velocity. Assuming i = 53° and n c ~ 10 -3 cm -2 
we derive from the [SII] A673.1nm flux a mass outflow rate 
of logM out [M Q /yr] ~ —10.4. Combining this with the con- 
temporaneous measurement of the mass accretion rate we find 
an outflow-to-inflow rate M out /M acc ~ 0.3. This value for 
the mass outflow rate of FU Tau A could be both a lower limit 
(given that n c is an upper limit) or an upper limit (given that 
i is a lower limit); see Sect. |4] for details. Very few measure- 
ments of Aiout /Macc have been presented in the literature for 
VLM objects. While some of these observations suggest that 
mass inflow and outflow rate are similar in the VLM regime 
dBacciotti et all 120111: IWhelan et all l2009h . others point at re- 
sults close to the 'canonical' value for higher-mass cTTS of 
M out /M acc ~ 0.1. ..0.01 dJoergens et al.ll2012l) . A final conclu- 
sion on the comparison of the mass outflow to inflow ratio be- 
tween BDs and cTTSs and the compari son of obs ervations to the 
predictions of jet-launching models (Cabrit 2009) is premature 
considering the large uncertainties associated with the observa- 
tions of M ut an d Ma.cc in VLM objects and the inhomogeneous 
approach by different authors for calculating these numbers. 

Signs for in- and outflows of FU Tau A are also evident from 
the line profiles and line shifts observed in the X-Shooter spec- 
trum. The modest deficiency of flux on the red side of the low- 
n Balmer lines is in qualitative agreement with the outcome of 
magnetosph eric accretion models for high inclination. The cal- 
culations bv lHartmann et all d 19941) predict a lack of flux in the 
red wing of optically thin emission lines that increases with de- 
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creasing inclination angle due to occultation of the receeding 
part of the accretion flow by the disk. The strongest forbidden 
lines show a similar blueshift of « — lOkm/s with respect to the 
stellar motion. This emission likely represents the approaching 
lobe of the jet, while the receeding (redshifted) part of the out- 
flow may be occulted by the accretion disk. These are plausible 
interpretations which are consistent with our range for the disk 
inclination angle (derived from v sin i, rotation period and stellar 
radius) that excludes a pole-on view. 

6.2. FU Tau A in the HR diagram 

6.2.1. Extinction 

A wrong extinction estimate could, in principle, misplace the 
object in the HR diagram in the vertical direction. We have 
determined the spectral type of FU Tau A to M6.5 ± 0.4 and 
the optical extinction to Ay = 0.5 ± .5 mag. This is slightl y 
hotter and less absorbed than found by Luhma net al.l (|2009a), 
and results in a higher effective temperature (T e g = 2940 K) 
and smaller bolometric luminosity (Lboi = 0.13 L ) and radius 
(i?* = 1.38 R©). The smaller extinction moves the object ver- 
tically down in the HR diagram by less than a factor two. This 
is by far not sufficient for resolving the luminosity problem of 
FUTauA. 

6.2.2. Accretion 

A possible explanation for the overluminosity of FU Tau A with 
respect to the youngest isochrones of evolutionary models could 
be an overestimate of the bolometric lumino sity due to a substan- 
tial contribution of accretion luminosity. In IScholz et al.1 (1201 2l) 
we have argued that this scenario is unlikely because of the rela- 
tively low temperature derived for the accretion hot spots. It was 
also pointed out that a significantly lower value for Lhn i , whe n 
combined with the X-ray data presented by Stel zer et al.l (12010), 
would yield an exceptionally high X-ray to bolometric flux ratio. 
The X-Shooter spectrum has now allowed us to directly compute 
the accretion luminosity by modeling of the Balmer continuum. 
We found that accretion makes up for only 0.1 % of the bolomet- 
ric luminosity of FU Tau A. Therefore, accretion luminosity can 
not explain the overluminosity of FU Tau A. This is also corrob- 
orated by the absence of strong veiling in the X-Shooter spec- 
trum. 

Accretion possibly plays a role in another puzzling prop- 
erty of FU Tau A, its weak lithium absorption (see discussion in 
Sect.lOTl. 

6.2.3. Magnetic activity 

An alternative possibility to bring FUTauA in better agree- 
ment with evolutionary models would be to shift the object 
horizontally to wards higher effectiv e te mperature and mass. It 
was sh own by IChabrier et al.l (120071) and iMacDonald & Mullanl 
(2009) that the influence of a strong magnetic field onto con- 
vection may provide such an effect. In this scenario FU Tau A 
appears cooler and with larger radius than expected for its mass. 
As T c g is obtained from observations without resorting to evolu- 
tionary models while mass is model-dependent, this would im- 
ply a higher mass for FU Tau A than the 0.08 M® obtained when 
ex trapolating the objec t vertically down to the 1 Myr isochrone 
of lBaraffeetandl998l). 

In Stelzer et al.1 d2010l) we had derived a mass of ~ 0.2 M 
for FU Tau A assuming this scenario based on the stellar param- 



eters given bv lLuhmanetail (|2009a) and we had estimated the 
rotation period and the magnetic field strength using standard 
magnetospheric accretion models. These estimates can now be 
updated making use of the new values for the stellar parameters 
derived from the X-Shooter spectrum. 

First, the new value for the bolometric luminosity corre- 
sponds to a mass of 0.1 5M^ and a te mperature of 3090 K on 
the 1 Myr isochrone of iBaraffe et al.l (119981) . Our new value 
for the mass accretion rate (see Sect. 15.2. I) was based on a 
mass of M — 0.08 M Q . This would be corrected downward 
to log (M acc ) [M /yr] = -10.1 ± 0.2. Using these new pa- 
rameters and an inner disk truncation r adius of R ln ~ 2. ..5 R* 
the magnetic field strength predicted by KoenigJ (Il99ll) is B ps 
40...180G. No magnetic field measurement has been performed 
yet for FUTauA that would allow to test these numbers. In any 
case, this estimate shows that the combination of stellar pa- 
rameters and mass accretion rate observed for FUTauA sup- 
port t he presence of a stron g surface field. The X-wind model 
(e.g.. iMohantv & Shul 120 08). where the Keplerian angular ve- 
locity is equal to the stellar rotational velocity at the inner disk 
truncation radius, predicts a stellar rotation period of P rot « 
1.4...5.5d for the range of R j n given above. The value ob- 
served by IScholz et al.l (120121) for the rotation period (P ro t = 
4.0 ± 0.2 d) is consistent with this prediction. As pointed out in 
our previous work this period is rather large for substellar mass 
according to the empiric al trend between period and mass (e.g., 
IScholz & EisiHffell200l . 

The influence of magnetic activity on the stellar parameters 
is also expected to have an effect on the lithium depletion of late- 
type stars but this is unlikely to be relevant for FU Tau A as we 
describe in Sect. 16.31 

6.2.4. Age 

Both FU Tau A and its brown dwarf companion FU Tau B have 
the highest absolute J band magnitude of all Taurus members 
with which they share the effectiv e temperature, T a ff ~ 2800 K 
and T off ~ 2400 K respectively dScholz et alj|2012|) . This sug- 
gests that the pair is coeval but younger than the other VLM ob- 
jects in Taurus. The low gravity of log g — 3.5 ± 0.5 measured 
from the X-Shooter spectrum is a confirmation of the pre-main 
sequence nature of FU Tau A. However, in the standard picture 
extreme youth is not compatible with the weak lithium signature. 
Furthermore, there are no signs for an envelope which would be 
expected if the object were still in a protostellar phase. 

6.2.5. Binarity 

The picture could be further complicated if FUTauA was a 
close binary. The maximum decrease in Lboi obtained in this 
scenario is a factor two which would make FU Tau A similar to 
some other V LM objects in Taurus that are overluminous with 
respect to the Baraffe et al. (1998) lMyr isochrone. An equal 
mass binary with 140 AU separation (« 1" at the distance of 
Taurus) would have an RV amplitude of less than ~ 1 km/s 
but a sub-arcsecond binary would produce a much larger RV 
signal. Binarity could be responsible for some of the difference 
between the RV we measured for FUTauA and most Taurus 
members. One of the most obvious effects of binarity, besides 
its influence on the RV, would be that the rotational velocity has 
an unknown contribution from the orbital motion and the mea- 
sured v sin i is an upper limit. The R sin i estimate derived from 
the rotational velocity and the rotation period is, in fact, slightly 
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larger than the stellar radius obtained from Stefan-Boltzmann's 
law(i?sinz = 1.6 ± 0.5 R©, i?* = 1.4 R Q ) but compatible with 
each other within the errors. The lithium equivalent width may 
also change in the case of binarity but in a way which is hard 
to predict without knowledge of the relative brightness of both 
binary components. 



6.2.6. Distance 

Finally, an obvious way to reduce the bolometric luminosity of 
FU Tau A would be a smaller distance. This would also lead to a 
better agreement of the X-ray luminosit y of FU Tau A with other 
objects of similar mass in Taurus (see IStelzer et al.lfeOlOb . For 
a smaller distance, the difference between the mean mass accre- 
tion rate obtained from the 24 emission lines and the mass accre- 
tion rate obtained from the Ha 10 % width of FU Tau A could be 
reconciled, e.g. a distance of 75 pc yields (log M acc ) [M Q /yr] = 
— 10.4 ± 0.2 for the mean of the emission lines with the stellar 
parameters from Sect. 13.21 The magnetic field would be larger 
(B = 90...430G) and the rotation faster (F rot = 0.8.. .2.9d) 
than predicted for the magnetic activity scenario of Sect. 16.2.31 
However, the theoretical predictions for these parameters should 
be taken as face values and not as firm prove for one or the other 
hypothesis. 

The closer distance poses the problem of either the existence 
of a dark cloud apparently associated with the Taurus filaments 
but in reality located in the foreground or the unlikely projec- 
tion of an isolated young binary onto a background dark cloud. 
An argument in favor of the latter case comes from the differ- 
ence of the LSR velocity of FU Tau A (12.5 ± 2.9 km/s) and the 
value derived for Barn ard 2 15 from CO maps (t>LSR ~ 7km/s; 
iNaravanan et al.ll2008l) . On the other hand, this scenario where 
FU Tau A is in the foreground of Barnard215, i.e. far from any 
signatures of star formation, puts a question mark on its origin. 
Ejection from a star forming site seems not a viable option as 
for typical ejection speeds of ~ lkm/s (e.g. iBate et al.l 12003) 
FU Tau A cannot have travelled be yond the (proje c ted) dia meter 
of Barnard215 within its lifetime. Luhma n et al. (2009a) men- 
tioned that the SDSS image shows nebulosity centered and sur- 
rounding FU Tau which, when combined with the signatures of 
youth, suggests that the object is associated with the dark cloud. 

The RV of FUTau A (22.5 ± 2.9 km/s) is marginally con- 
sistent with the RV distribution of Taurus and so is its UVW 
space velocity (({/, V, W) = (-22.09,-9 . 08,- 10 .36)km/s). For 
the proper motion from iLuhman et al.l (12009 al) a closer dis- 
tance results in a space motion with a stronger deviation 
from the mean Taurus UVW velocities (e.g., (U, V, W) = 
(—21.80,— 3.98,— 8. 64km/s) for 80 pc). However, considering 
the uncertainty of the proper motion inferred from the different 
literature sources, the change in space motion associated with a 
closer distance seems not significant enough to allow us to draw 
conclusions. 



6.3. The lithium problem 

The lithium equivalent width measured in the X-Shooter spec- 
trum of FU Tau A is ~ 430 mA. This corresponds to a lithium 
abundance below A (Li) ~ 1.9 according to th e curves of 
growth calculated by Zapate ro Osorio et al. (120021) . At present, 
FUTau A defines the lower envelope in the lithium equiva- 
lent widths of spectral types M6...M7. We have shown that 
this does not change even if weak veiling, consistent with the 
depth of the TiO absorption bands, is assumed. The second low- 



est value is represented by KPNO-T au5, another young BD 
in Tau rus (Bar radq v Navascu es 2004). This object has no disk 
dMonin et alj|2010l) and is below the canoni cal threshold for ac- 
cretor s in terms of its Ha equivalent width dBarrado v Navascuesl 
2004) such that veiling seems unlikely to be present. Lithium 
measurements are available for a limited sample of such late- 
type young objects. Therefore, no conclusions can be drawn yet 
on the typical value at those spectral types. However, the large 
dispersion of lithium abundances is seen also at early-M types 
and has been noted in the literature as an unsolved problem. 

Quite a number of parameters not taken into account in stan- 
dard PMS ev olutionary models may af fect the lithium deple- 
tion timescale. Egg enberger et al.l (120121) described the influence 
of rotation and disk lifetime on the lithium depletion of solar- 
mass stars. Interestingly, rotation shifts the evolutionary tracks 
to slightly lower T e g in the HR diagram. Similarly, if the disks 
are longer lived lithium gets more and more rapidly depleted. 
The impact of rotation and disk lifetime on lithium depletion is 
at odds with the fact that, for a given star, ages derived from its 
lithium abundance tend to be older than those obt ained from its 
position in the HR diagram. lYee & Jensen! d2010t) explain how 
the HR diagram ages and the lithium ages can be brought into 
qualitative agreement by the inflation of the radii due to mag- 
netic activity described in Sect. 16.2.31 The cooler effective tem- 
peratures associated with the larger radius would imply for a 
given star an older age from the HR diagram and a younger 
lithium depletion age. However, all these studies regard solar- 
type stars and their validity in the VLM regime can not be taken 
for granted. Moreover, these arguments should break down for 
objects as young as FU Tau A where no lithium depletion at all 
is expected from pre-main sequence models. 

For FU Tau A anothe r scenario is more plausible. 
Bara ffe & Chabrierl (1201 Ob argue that episodic accretion 
bursts during the pre-main sequence phase strongly decrease the 
timescale for lithium depletion in low-mass stars. In particular, 
they show that an accreting 0.1 M Q star may completely deplete 
its initial lithium content within 10 Myrs compared to > 50 Myrs 
required in the non-accreting case. For suitable values of the 
initial protostar mass and the strength and number of accretion 
outbursts the observed lithium abundance of FU Tau A seems 
in agreement with the predicted abundance from these models . 
However, in strong contrast to the observation for FUTau A, 
the inclusion of accretion in evolutionar y models produce s 
an underluminosity in the HR diagram dBaraffe et al.l l2009h . 
Special cases of these models evolve through a high-luminosity 
phase but this phase is very short-lived and the probability 
of observing its realization in nature is rather low. Moreover, 
FUTau A is the only known VLM YSO with strong lithium 
under-abundance that is a confirmed accretor. 

A strong underabundance of lithiu m with respec t to the cos- 
mic abundance has also been found bv lJohnas et al.l d200jj) who 
modeled the high-resolution spectra of several young VLM ob- 
jects in the Cha I star forming region with synthetic atmospheres. 
Their targets have spectral types M6/M7 and lithium abundances 
A(Li) < 2.0, just as FUTau A. This seems to imply that our un- 
derstanding of lithium depletion or line formation in VLM ob- 
jects is incomplete. 

7. Summary and conclusions 

We have analysed in detail an X-Shooter spectrum of FU Tau A 
which has allowed us to redetermine its stellar parameters and 
has provided a wealth of new information on activity, accre- 
tion and outflow diagnostics. This allows us to exclude accre- 
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tion luminosity and extinction as cause of the overluminosity of 
FU Tau A in the HR diagram. An inflated radius as a result of 
magnetic activity leading to an underestimate of T e g and M„ 
might be responsible for it. Strong activity is consistent with the 
high X-ray luminosity of FU Tau A. A combination of magnetic 
activity and very young age for the FU Tau binary may be re- 
sponsible for its HR diagram location. Another possibility is a 
closer distance but this is hard to explain due to the association 
of FU Tau A with the Barnard 215 dark cloud that is believed to 
be part of the Taurus complex. The possibility of FU Tau A be- 
ing an unresolved binary must be taken into account. This would 
bring FU Tau A in better agreement with other luminous Taurus 
members and possibly with the RV distribution of Taurus. High- 
resolution spectroscopy is needed to search for a putative com- 
panion and to verify the unexpected evidence for lithium deple- 
tion. However, even with the available information the startling 
properties of FU Tau A makes it a benchmark object to push on 
the further development of PMS evolutionary models. 
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